\(\int \frac {\cos ^2(e+f x) (a+a \sin (e+f x))^{7/2}}{(c-c \sin (e+f x))^{17/2}} \, dx\) [43]

Optimal result

Integrand size = 38, antiderivative size = 145 \[ \int \frac {\cos ^2(e+f x) (a+a \sin (e+f x))^{7/2}}{(c-c \sin (e+f x))^{17/2}} \, dx=\frac {\cos (e+f x) (a+a \sin (e+f x))^{9/2}}{14 a c f (c-c \sin (e+f x))^{15/2}}+\frac {\cos (e+f x) (a+a \sin (e+f x))^{9/2}}{84 a c^2 f (c-c \sin (e+f x))^{13/2}}+\frac {\cos (e+f x) (a+a \sin (e+f x))^{9/2}}{840 a c^3 f (c-c \sin (e+f x))^{11/2}} \]

[Out]

1/14*cos(f*x+e)*(a+a*sin(f*x+e))^(9/2)/a/c/f/(c-c*sin(f*x+e))^(15/2)+1/84*cos(f*x+e)*(a+a*sin(f*x+e))^(9/2)/a/
c^2/f/(c-c*sin(f*x+e))^(13/2)+1/840*cos(f*x+e)*(a+a*sin(f*x+e))^(9/2)/a/c^3/f/(c-c*sin(f*x+e))^(11/2)

Rubi [A] (verified)

Time = 0.42 (sec) , antiderivative size = 145, normalized size of antiderivative = 1.00, number of steps used = 4, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.079, Rules used = {2920, 2822, 2821} \[ \int \frac {\cos ^2(e+f x) (a+a \sin (e+f x))^{7/2}}{(c-c \sin (e+f x))^{17/2}} \, dx=\frac {\cos (e+f x) (a \sin (e+f x)+a)^{9/2}}{840 a c^3 f (c-c \sin (e+f x))^{11/2}}+\frac {\cos (e+f x) (a \sin (e+f x)+a)^{9/2}}{84 a c^2 f (c-c \sin (e+f x))^{13/2}}+\frac {\cos (e+f x) (a \sin (e+f x)+a)^{9/2}}{14 a c f (c-c \sin (e+f x))^{15/2}} \]

[In]

Int[(Cos[e + f*x]^2*(a + a*Sin[e + f*x])^(7/2))/(c - c*Sin[e + f*x])^(17/2),x]

[Out]

(Cos[e + f*x]*(a + a*Sin[e + f*x])^(9/2))/(14*a*c*f*(c - c*Sin[e + f*x])^(15/2)) + (Cos[e + f*x]*(a + a*Sin[e
+ f*x])^(9/2))/(84*a*c^2*f*(c - c*Sin[e + f*x])^(13/2)) + (Cos[e + f*x]*(a + a*Sin[e + f*x])^(9/2))/(840*a*c^3
*f*(c - c*Sin[e + f*x])^(11/2))

Rule 2821

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Simp
[b*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*((c + d*Sin[e + f*x])^n/(a*f*(2*m + 1))), x] /; FreeQ[{a, b, c, d, e, f
, m, n}, x] && EqQ[b*c + a*d, 0] && EqQ[a^2 - b^2, 0] && EqQ[m + n + 1, 0] && NeQ[m, -2^(-1)]

Rule 2822

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Simp
[b*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*((c + d*Sin[e + f*x])^n/(a*f*(2*m + 1))), x] + Dist[(m + n + 1)/(a*(2*m
+ 1)), Int[(a + b*Sin[e + f*x])^(m + 1)*(c + d*Sin[e + f*x])^n, x], x] /; FreeQ[{a, b, c, d, e, f, m, n}, x]
&& EqQ[b*c + a*d, 0] && EqQ[a^2 - b^2, 0] && ILtQ[Simplify[m + n + 1], 0] && NeQ[m, -2^(-1)] && (SumSimplerQ[m
, 1] ||  !SumSimplerQ[n, 1])

Rule 2920

Int[cos[(e_.) + (f_.)*(x_)]^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*((c_) + (d_.)*sin[(e_.) + (f_.)*
(x_)])^(n_.), x_Symbol] :> Dist[1/(a^(p/2)*c^(p/2)), Int[(a + b*Sin[e + f*x])^(m + p/2)*(c + d*Sin[e + f*x])^(
n + p/2), x], x] /; FreeQ[{a, b, c, d, e, f, n, p}, x] && EqQ[b*c + a*d, 0] && EqQ[a^2 - b^2, 0] && IntegerQ[p
/2]

Rubi steps \begin{align*} \text {integral}& = \frac {\int \frac {(a+a \sin (e+f x))^{9/2}}{(c-c \sin (e+f x))^{15/2}} \, dx}{a c} \\ & = \frac {\cos (e+f x) (a+a \sin (e+f x))^{9/2}}{14 a c f (c-c \sin (e+f x))^{15/2}}+\frac {\int \frac {(a+a \sin (e+f x))^{9/2}}{(c-c \sin (e+f x))^{13/2}} \, dx}{7 a c^2} \\ & = \frac {\cos (e+f x) (a+a \sin (e+f x))^{9/2}}{14 a c f (c-c \sin (e+f x))^{15/2}}+\frac {\cos (e+f x) (a+a \sin (e+f x))^{9/2}}{84 a c^2 f (c-c \sin (e+f x))^{13/2}}+\frac {\int \frac {(a+a \sin (e+f x))^{9/2}}{(c-c \sin (e+f x))^{11/2}} \, dx}{84 a c^3} \\ & = \frac {\cos (e+f x) (a+a \sin (e+f x))^{9/2}}{14 a c f (c-c \sin (e+f x))^{15/2}}+\frac {\cos (e+f x) (a+a \sin (e+f x))^{9/2}}{84 a c^2 f (c-c \sin (e+f x))^{13/2}}+\frac {\cos (e+f x) (a+a \sin (e+f x))^{9/2}}{840 a c^3 f (c-c \sin (e+f x))^{11/2}} \\ \end{align*}

Mathematica [B] (verified)

Leaf count is larger than twice the leaf count of optimal. \(419\) vs. \(2(145)=290\).

Time = 12.87 (sec) , antiderivative size = 419, normalized size of antiderivative = 2.89 \[ \int \frac {\cos ^2(e+f x) (a+a \sin (e+f x))^{7/2}}{(c-c \sin (e+f x))^{17/2}} \, dx=\frac {16 \left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right )^3 (a (1+\sin (e+f x)))^{7/2}}{7 f \left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right )^7 (c-c \sin (e+f x))^{17/2}}-\frac {16 \left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right )^5 (a (1+\sin (e+f x)))^{7/2}}{3 f \left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right )^7 (c-c \sin (e+f x))^{17/2}}+\frac {24 \left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right )^7 (a (1+\sin (e+f x)))^{7/2}}{5 f \left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right )^7 (c-c \sin (e+f x))^{17/2}}-\frac {2 \left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right )^9 (a (1+\sin (e+f x)))^{7/2}}{f \left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right )^7 (c-c \sin (e+f x))^{17/2}}+\frac {\left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right )^{11} (a (1+\sin (e+f x)))^{7/2}}{3 f \left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right )^7 (c-c \sin (e+f x))^{17/2}} \]

[In]

Integrate[(Cos[e + f*x]^2*(a + a*Sin[e + f*x])^(7/2))/(c - c*Sin[e + f*x])^(17/2),x]

[Out]

(16*(Cos[(e + f*x)/2] - Sin[(e + f*x)/2])^3*(a*(1 + Sin[e + f*x]))^(7/2))/(7*f*(Cos[(e + f*x)/2] + Sin[(e + f*
x)/2])^7*(c - c*Sin[e + f*x])^(17/2)) - (16*(Cos[(e + f*x)/2] - Sin[(e + f*x)/2])^5*(a*(1 + Sin[e + f*x]))^(7/
2))/(3*f*(Cos[(e + f*x)/2] + Sin[(e + f*x)/2])^7*(c - c*Sin[e + f*x])^(17/2)) + (24*(Cos[(e + f*x)/2] - Sin[(e
+ f*x)/2])^7*(a*(1 + Sin[e + f*x]))^(7/2))/(5*f*(Cos[(e + f*x)/2] + Sin[(e + f*x)/2])^7*(c - c*Sin[e + f*x])^
(17/2)) - (2*(Cos[(e + f*x)/2] - Sin[(e + f*x)/2])^9*(a*(1 + Sin[e + f*x]))^(7/2))/(f*(Cos[(e + f*x)/2] + Sin[
(e + f*x)/2])^7*(c - c*Sin[e + f*x])^(17/2)) + ((Cos[(e + f*x)/2] - Sin[(e + f*x)/2])^11*(a*(1 + Sin[e + f*x])
)^(7/2))/(3*f*(Cos[(e + f*x)/2] + Sin[(e + f*x)/2])^7*(c - c*Sin[e + f*x])^(17/2))

Maple [A] (verified)

Time = 0.39 (sec) , antiderivative size = 200, normalized size of antiderivative = 1.38

method	result	size
default	\(\frac {\sqrt {a \left (1+\sin \left (f x +e \right )\right )}\, a^{3} \left (9 \left (\cos ^{5}\left (f x +e \right )\right ) \sin \left (f x +e \right )-63 \left (\cos ^{5}\left (f x +e \right )\right )-216 \left (\cos ^{3}\left (f x +e \right )\right ) \sin \left (f x +e \right )+469 \left (\cos ^{3}\left (f x +e \right )\right )+790 \cos \left (f x +e \right ) \sin \left (f x +e \right )-854 \cos \left (f x +e \right )-688 \tan \left (f x +e \right )+448 \sec \left (f x +e \right )\right )}{105 f \left (\cos ^{6}\left (f x +e \right )+6 \left (\cos ^{4}\left (f x +e \right )\right ) \sin \left (f x +e \right )-18 \left (\cos ^{4}\left (f x +e \right )\right )-32 \left (\cos ^{2}\left (f x +e \right )\right ) \sin \left (f x +e \right )+48 \left (\cos ^{2}\left (f x +e \right )\right )+32 \sin \left (f x +e \right )-32\right ) \sqrt {-c \left (\sin \left (f x +e \right )-1\right )}\, c^{8}}\)	\(200\)

[In]

int(cos(f*x+e)^2*(a+a*sin(f*x+e))^(7/2)/(c-c*sin(f*x+e))^(17/2),x,method=_RETURNVERBOSE)

[Out]

1/105/f*(a*(1+sin(f*x+e)))^(1/2)*a^3/(cos(f*x+e)^6+6*cos(f*x+e)^4*sin(f*x+e)-18*cos(f*x+e)^4-32*cos(f*x+e)^2*s
in(f*x+e)+48*cos(f*x+e)^2+32*sin(f*x+e)-32)/(-c*(sin(f*x+e)-1))^(1/2)/c^8*(9*cos(f*x+e)^5*sin(f*x+e)-63*cos(f*
x+e)^5-216*cos(f*x+e)^3*sin(f*x+e)+469*cos(f*x+e)^3+790*cos(f*x+e)*sin(f*x+e)-854*cos(f*x+e)-688*tan(f*x+e)+44
8*sec(f*x+e))

Fricas [A] (verification not implemented)

none

Time = 0.34 (sec) , antiderivative size = 205, normalized size of antiderivative = 1.41 \[ \int \frac {\cos ^2(e+f x) (a+a \sin (e+f x))^{7/2}}{(c-c \sin (e+f x))^{17/2}} \, dx=-\frac {{\left (35 \, a^{3} \cos \left (f x + e\right )^{4} - 154 \, a^{3} \cos \left (f x + e\right )^{2} + 128 \, a^{3} - 14 \, {\left (5 \, a^{3} \cos \left (f x + e\right )^{2} - 8 \, a^{3}\right )} \sin \left (f x + e\right )\right )} \sqrt {a \sin \left (f x + e\right ) + a} \sqrt {-c \sin \left (f x + e\right ) + c}}{105 \, {\left (7 \, c^{9} f \cos \left (f x + e\right )^{7} - 56 \, c^{9} f \cos \left (f x + e\right )^{5} + 112 \, c^{9} f \cos \left (f x + e\right )^{3} - 64 \, c^{9} f \cos \left (f x + e\right ) - {\left (c^{9} f \cos \left (f x + e\right )^{7} - 24 \, c^{9} f \cos \left (f x + e\right )^{5} + 80 \, c^{9} f \cos \left (f x + e\right )^{3} - 64 \, c^{9} f \cos \left (f x + e\right )\right )} \sin \left (f x + e\right )\right )}} \]

[In]

integrate(cos(f*x+e)^2*(a+a*sin(f*x+e))^(7/2)/(c-c*sin(f*x+e))^(17/2),x, algorithm="fricas")

[Out]

-1/105*(35*a^3*cos(f*x + e)^4 - 154*a^3*cos(f*x + e)^2 + 128*a^3 - 14*(5*a^3*cos(f*x + e)^2 - 8*a^3)*sin(f*x +
e))*sqrt(a*sin(f*x + e) + a)*sqrt(-c*sin(f*x + e) + c)/(7*c^9*f*cos(f*x + e)^7 - 56*c^9*f*cos(f*x + e)^5 + 11
2*c^9*f*cos(f*x + e)^3 - 64*c^9*f*cos(f*x + e) - (c^9*f*cos(f*x + e)^7 - 24*c^9*f*cos(f*x + e)^5 + 80*c^9*f*co
s(f*x + e)^3 - 64*c^9*f*cos(f*x + e))*sin(f*x + e))

Sympy [F(-1)]

Timed out. \[ \int \frac {\cos ^2(e+f x) (a+a \sin (e+f x))^{7/2}}{(c-c \sin (e+f x))^{17/2}} \, dx=\text {Timed out} \]

[In]

integrate(cos(f*x+e)**2*(a+a*sin(f*x+e))**(7/2)/(c-c*sin(f*x+e))**(17/2),x)

[Out]

Timed out

Maxima [F(-1)]

Timed out. \[ \int \frac {\cos ^2(e+f x) (a+a \sin (e+f x))^{7/2}}{(c-c \sin (e+f x))^{17/2}} \, dx=\text {Timed out} \]

[In]

integrate(cos(f*x+e)^2*(a+a*sin(f*x+e))^(7/2)/(c-c*sin(f*x+e))^(17/2),x, algorithm="maxima")

[Out]

Timed out

Giac [A] (verification not implemented)

none

Time = 0.32 (sec) , antiderivative size = 202, normalized size of antiderivative = 1.39 \[ \int \frac {\cos ^2(e+f x) (a+a \sin (e+f x))^{7/2}}{(c-c \sin (e+f x))^{17/2}} \, dx=-\frac {{\left (35 \, a^{3} \sqrt {c} \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{8} - 105 \, a^{3} \sqrt {c} \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{6} + 126 \, a^{3} \sqrt {c} \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{4} - 70 \, a^{3} \sqrt {c} \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} + 15 \, a^{3} \sqrt {c} \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right )\right )} \sqrt {a}}{840 \, c^{9} f \mathrm {sgn}\left (\sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) \sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{14}} \]

[In]

integrate(cos(f*x+e)^2*(a+a*sin(f*x+e))^(7/2)/(c-c*sin(f*x+e))^(17/2),x, algorithm="giac")

[Out]

-1/840*(35*a^3*sqrt(c)*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e))*sin(-1/4*pi + 1/2*f*x + 1/2*e)^8 - 105*a^3*sqrt(c)*
sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e))*sin(-1/4*pi + 1/2*f*x + 1/2*e)^6 + 126*a^3*sqrt(c)*sgn(cos(-1/4*pi + 1/2*f
*x + 1/2*e))*sin(-1/4*pi + 1/2*f*x + 1/2*e)^4 - 70*a^3*sqrt(c)*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e))*sin(-1/4*pi
+ 1/2*f*x + 1/2*e)^2 + 15*a^3*sqrt(c)*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e)))*sqrt(a)/(c^9*f*sgn(sin(-1/4*pi + 1
/2*f*x + 1/2*e))*sin(-1/4*pi + 1/2*f*x + 1/2*e)^14)

Mupad [B] (verification not implemented)

Time = 15.98 (sec) , antiderivative size = 764, normalized size of antiderivative = 5.27 \[ \int \frac {\cos ^2(e+f x) (a+a \sin (e+f x))^{7/2}}{(c-c \sin (e+f x))^{17/2}} \, dx=\text {Too large to display} \]

[In]

int((cos(e + f*x)^2*(a + a*sin(e + f*x))^(7/2))/(c - c*sin(e + f*x))^(17/2),x)

[Out]

-((c - c*((exp(- e*1i - f*x*1i)*1i)/2 - (exp(e*1i + f*x*1i)*1i)/2))^(1/2)*((a^3*exp(e*4i + f*x*4i)*(a + a*((ex
p(- e*1i - f*x*1i)*1i)/2 - (exp(e*1i + f*x*1i)*1i)/2))^(1/2)*16i)/(3*c^9*f) + (64*a^3*exp(e*5i + f*x*5i)*(a +
a*((exp(- e*1i - f*x*1i)*1i)/2 - (exp(e*1i + f*x*1i)*1i)/2))^(1/2))/(3*c^9*f) - (a^3*exp(e*6i + f*x*6i)*(a + a
*((exp(- e*1i - f*x*1i)*1i)/2 - (exp(e*1i + f*x*1i)*1i)/2))^(1/2)*1088i)/(15*c^9*f) - (576*a^3*exp(e*7i + f*x*
7i)*(a + a*((exp(- e*1i - f*x*1i)*1i)/2 - (exp(e*1i + f*x*1i)*1i)/2))^(1/2))/(5*c^9*f) + (a^3*exp(e*8i + f*x*8
i)*(a + a*((exp(- e*1i - f*x*1i)*1i)/2 - (exp(e*1i + f*x*1i)*1i)/2))^(1/2)*5472i)/(35*c^9*f) + (576*a^3*exp(e*
9i + f*x*9i)*(a + a*((exp(- e*1i - f*x*1i)*1i)/2 - (exp(e*1i + f*x*1i)*1i)/2))^(1/2))/(5*c^9*f) - (a^3*exp(e*1
0i + f*x*10i)*(a + a*((exp(- e*1i - f*x*1i)*1i)/2 - (exp(e*1i + f*x*1i)*1i)/2))^(1/2)*1088i)/(15*c^9*f) - (64*
a^3*exp(e*11i + f*x*11i)*(a + a*((exp(- e*1i - f*x*1i)*1i)/2 - (exp(e*1i + f*x*1i)*1i)/2))^(1/2))/(3*c^9*f) +
(a^3*exp(e*12i + f*x*12i)*(a + a*((exp(- e*1i - f*x*1i)*1i)/2 - (exp(e*1i + f*x*1i)*1i)/2))^(1/2)*16i)/(3*c^9*
f)))/(exp(e*1i + f*x*1i)*14i - 90*exp(e*2i + f*x*2i) - exp(e*3i + f*x*3i)*350i + 910*exp(e*4i + f*x*4i) + exp(
e*5i + f*x*5i)*1638i - 2002*exp(e*6i + f*x*6i) - exp(e*7i + f*x*7i)*1430i - exp(e*9i + f*x*9i)*1430i + 2002*ex
p(e*10i + f*x*10i) + exp(e*11i + f*x*11i)*1638i - 910*exp(e*12i + f*x*12i) - exp(e*13i + f*x*13i)*350i + 90*ex
p(e*14i + f*x*14i) + exp(e*15i + f*x*15i)*14i - exp(e*16i + f*x*16i) + 1)